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Faint K-Selected Galaxy Correlations and Clustering Evolution 

R. G. Carlberg, 1 Lennox L. Cowie, 2,3 ' 4 Antoinette Songaila, 2 ' 3 and Esther M. Hu 2 ' 3 ' 4 

ABSTRACT 

Angular and spatial correlations are measured for i^-band-selected galaxies, 248 
having redshifts, 54 with z > 1, in two patches of combined area ~ 27 arcmin 2 . The 
angular correlation for K < 21. 5 mag is w(0) ~ (9/1A ± 0.19" e ±tu )-°- 8 . From the 
redshift sample we find that the real-space correlation, calculated with qo = 0.1, of 
Mk < —23.5 mag galaxies (k-corrected) is £(r) = (r/2.9e ±012 /i~ 1 Mpc) -1 ' 8 at a mean 
z ~ 0.34, (r/2.0e ±0 - 15 /i- 1 Mpc)- L8 at z ~ 0.62, (r/l^e^- 15 ^ 1 Mpc) -1 - 8 at z ~ 0.97, 
and (r /l.Oe^^h -1 Mpc) -1 ' 8 at z ~ 1.39, the last being a formal upper limit for a blue- 
biased sample. In general, these are more correlated than optically selected samples in 
the same redshift ranges. Over the interval 0.3 < z < 0.9 galaxies with red rest-frame 
colors, (U — K)q > 2 AB mag, have £(r) ~ (r/2.4e ±ai4 /i -1 Mpc)~ L8 whereas bluer 
galaxies, which have a mean B of 23.7 mag and mean [O n] equivalent width W eq = 41 A, 
are very weakly correlated, with £(r) ~ (r/0.9e ±0 ' 22 /i _1 Mpc) -1 ' 8 . For our measured 
growth rate of clustering, this blue population, if non-merging, can grow only into a 
low-redshift population less luminous than 0.4L*. The cross-correlation of low- and high- 
luminosity galaxies at z ~ 0.6 appears to have an excess in the correlation amplitude 
within 100 /i -1 kpc. The slow redshift evolution is consistent with these galaxies tracing 
the mass clustering in low density, O ~ 0.2, relatively unbiased, ag — 0.8, universe, but 
cannot yet exclude other possibilities. 
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1. Introduction 



N-body simulations give reliable predictions for the redshift dependence of the two-point cor- 
relation function of the density field, £(r\z), as a function of f2. A convenient power-law parame- 



terization to describe the evolving correlation function of galaxies is ( Groth fc Peebles 1977 , Koo 



fc Szalay 1981 



C(r\z) = - 



r 



:i+z) 



-(3+6) 



(1) 



where the lengths r are measured in physical (proper) co-ordinates. For this double power-law 
approximation the predicted evolution of clustering in the mass field is faster for £1 = 1 (e = 1.0±0.1) 
than it is for low-f2 values, for instance e = 0.2 ± 0.1 for f2 = 0.2 (Colin, Carlberg & Couchman 
1996). Therefore, measurement of the redshift evolution of clustering can be used to test the 



gravitational instability theory of structure formation and the relation of galaxy clustering to dark 
matter clustering, and provides a constraint on £1. Knowledge of these quantities enables predictions 
of the distribution of assembly times of dark halos, which on the relatively small scales investigated 
here is of great interest for the mass evolution of galaxies. 

At present, observational measures of clustering evolution are uncertain simply due to the 
difficulties of assembling large samples of faint galaxies with sufficient sky coverage to give a statis- 
tically representative sample. At low redshift the form of nonlinear galaxy clustering is accurately 
established (e.g. Davis Sz Peebles 1983 , Loveday et al. 1995| , [Lin 1996] , Tucker et al. 1996 ), 



with a basic characterization of its dependencies on galaxy color, luminosity, and morphology. At 
higher redshifts the clustering is only now being directly measured ( Le Fevre et al. 1996j , [Shepherd 



et al. 1997), although the small fields leave concerns that field-to-field variations are not yet well 
controlled. 

The galaxy luminosity function and its color dependence evolve substantially over the redshift 
to 1 interval ( |Lilly et al. 1995| , [Ellis et al. 1~996| , |Cowie et al. 1996[ [Lin et al. 1997[ ). Differential 
luminosity evolution of blue and red galaxies, in which the blue galaxies are less correlated at 
low redshift, can cause the apparent correlation of a magnitude-limited sample to change faster 
than either of the two underlying populations are changing. In this paper we report the clustering 
properties of a very deep redshift survey selected in the K band. A near-IR selected survey has the 
enormous advantage that both k-corrections and the evolutionary corrections are small, allowing 
galaxy luminosities to be identified with total stellar mass with reasonable confidence. The Hawaii 
-ftT-band survey ( powie et al. 1996| ) with a couple hundred galaxy redshifts, is large enough to be 
useful for correlation studies. Furthermore this survey contains galaxies up to a redshift of 2.19, 
which provides a fairly large redshift baseline over which correlation changes can be measured. 

The next section summarizes the sample properties. Measures of the angular correlation are 
given in Section 3 and of the real space correlation function in Section 4. In Section 5 the correlations 
of red and blue galaxies and of low- and high-luminosity galaxies are compared. Section 6 discusses 
the redshift evolution of galaxy correlations and compares the available data to various model 
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predictions. Section 7 summarizes our results. All measurements in this paper assume Hq = 
100 h^ 1 kms -1 Mpc -1 and qo = 0.1. It should be noted that correlation amplitudes at z — 1 are 
reduced by about 30% for q = 0.5. 



2. The Hawaii iT-Band Sample 

The Hawaii iT-selected redshift survey of two fields constitutes a nearly complete sample down 
to K = 20, I = 23, and B = 24.5 mag. The survey is described in detail elsewhere (Cowie 



et al. 1996) although this analysis uses 20 new redshifts that have been recently obtained to 
complete the B selected subsample. The sample of 248 galaxies with redshifts in an area of about 
27 square arcminutes compares favourably with the moderate redshift CNOC sample, about 200 



galaxies in a single field covering 221 square arcminutes flShepherd et al. 1997 ), and with the 



591 galaxies covering about 71 square arcminutes (a result of the high sampling rate in 0.5' x 9.4' 
strips) in the CFRS study ( LeFevre et al. 1996] ) . On the other hand, the redshifts here extend 



all the way from 0.08 to 2.19, with 80% between 0.28 and 1.39. The number in any one redshift 
interval is sufficient to make useful correlation measurements (Figure |l]). All magnitudes and colors 
used in this paper are k-corrected. The sky positions of the objects are plotted in Figure ^, where 
open squares are galaxies, crosses represent stars, and objects that were not observed or which lack 
secure redshift identifications are shown with triangles. 

The galaxy redshifts, plotted against the projected physical distance in the RA direction from 
the field center, are shown in Figure |3|. The symbol area is proportional to the luminosity of the 
galaxy. The sample is known to be somewhat incomplete for the redder galaxies at the faintest 
magnitudes, which are generally expected to be galaxies beyond redshift one with low star formation 
rates. Incomplete samples, provided that they have no spatial bias, do not pose a problem for 
correlation studies, provided that the unclustered background distribution is generated from a 
smoothed version of the observed redshift distribution. We will approximate the observed redshift 
distributions as being constant over the various subranges of interest, which tests indicate to be an 
adequate approximation for these data. 



3. The Angular Correlation Function 

The angular correlation function is estimated as uj(6) = (DD — 2DR + RR)/RR ( [Landy fc 
Szalay 1993| ) where in a given range of angles DD is the number of data pairs, DR is the number of 
pairs between the data sample and a uniform random sample. This estimator is particularly useful 
when the clustering amplitude is significantly less than unity, as it is here. The resulting angular 
correlation function for the full photometric sample, K < 21.5 mag, is shown in Figure The 
errors are assigned using the bootstrap method ( [Efron Tibshirani 1986 ) with 100 resamplings. 
The angular galaxy correlations are diluted by the uncorrelated foreground stars, because the full 
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photometric sample is used, which requires that the correlation amplitude be corrected upward by 
a factor of (1 — f*)~ 2 , where /* is the fraction of objects which are stars in the two combined fields. 
We estimate that = 0.25, based on the stellar fraction of the spectroscopically identified sample 
that are stars. The uncorrected correlation function is lo(6) ~ (#/0"75)~ ' 8 , from which we find 
that the correlation angle of the galaxies is 0q = l."4±0."2 arcseconds, where the error is the internal 
error of the fit. Several recent studies in the optical region have found correlation angles of 8q — 1" 
at R = 23.5 mag flHudon fc Lilly 1996f) , O ^ 0"3 at I = 22.5 mag QLidman fc Peterson 1996[) and 



in a much fainter sample to r = 26 mag, 8q ~ 0"06 ( Brainerd, Smail &: Mould 1995 ). Since our 
mean (I — K) ~ 2.5 we note that the angular correlation that we observe is substantially larger 
than for a comparable I- or i?-band sample. This is a consequence of the K-band sample's being 
fundamentally more correlated than optical samples; we provide direct evidence for this below. 

In a survey done with a multi-object spectrograph there is the possibility that the angular 
correlation of the galaxies with redshifts is biased relative to the parent sample as a result of 
instrumental constraints for the selection of objects for spectroscopy. The simplest test for a 
selection effect of this type is to measure the ratio of DD pairs as a function of separation in the 
photometric and redshift sample as normalized to the total numbers in the two samples. We find 
that there is a 10% reduction in the number of pairs within about 10", and a 5% overselection 
at separations around 30", diminishing with increasing angle. This bias is generally less than the 
statistical errors, so we do not apply any geometric corrections in the correlation measurements. 



4. Real Space Correlations 

The velocity precision in this redshift survey is not adequate to allow a measurement of the 
redshift space correlation function at small scales. However, we can measure the projected real 
space correlation function, 



w p {Tv)= ZUr$ + y 2 )dy, (2) 

J — OO 

(Peebles 1980, Davis & Peebles 1983). The primary drawback of this estimator is that it averages 



over long cylindrical shells in redshift space, so it has a very broad window function ( peacock 1996 ) 



The real space correlation is derived from a power-law fit to the projected correlation function, 
Eq. HI, calculated from the galaxy sample. Operationally, w p (r p ) is the integral over r v of the 2D 
correlation function £(r p ,r v ), where r p and r v are the proper separation of galaxy pairs in the 
projected and redshift directions, respectively. Because our sample is not very large the errors are 
dominated by the small number statistics, rather than any complications of the estimator. We 
estimate £(r p ,r v ) as DD/ DR — 1, in the (r p ,r v ) co-ordinates. The random sample is 5 to 10xl0 4 
points. The sum over the r v co-ordinate is cut off at a practical range of 10/i -1 Mpc in proper 



co-ordinates at the redshift of the object (Shepherd et al. 1997). This value was selected from a 



range of trial values as being approximately the optimal value to maximize the signal-to-noise. This 
cutoff distance is sufficiently large that the correction to the integral for the correlation beyond the 
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cutoff is everywhere less than about 20%, and is ignored. The redshift distribution of the random 
sample is generated assuming that the unclustered redshift distribution is constant over the various 
redshift and luminosity ranges. Altering the redshift ranges shows this to be an entirely adequate 
approximation for these data. The errors are bootstrap estimates, reduced by a factor of \/3 (Mo 
Jing k Borner 1992|). 



The w p (r p ) measured over the redshift ranges 0.2-0.4, 0.4-0.8, 0.8-1.2, and 1.2-1.6 are shown 
in Figure |5[ The galaxies have minimum luminosities of Mk < —23.5 mag in the three higher 
redshift bins and Mk < —21.5 in the lower redshift range, to increase the sample size. A rest-frame 
luminosity of Mk = —23.5 mag is about half of L*. 



The w p (r p ) are fitted to a power-law correlation function, 

r(l)r(3-±) 



( Peebles 1980| ). The Gamma function factor is 3.68 for 7 = 1.8. Over the redshift range 0.2 < 



z < 0.4 we find £(r) ~ (r/2.9e ±ai2 /i~ 1 Mpc)" 1 ' 8 at a mean z ~ 0.34. For the 0.4 < z < 0.8 
range we find £(r) ~ (r/2.0e ±0 ' 15 /i _1 Mpc)~ L8 at a mean z ~ 0.62, and for 0.8 < z < 1.2, £(r) ~ 
(r/l^ - 15 ^- 1 Mpc) -1 - 8 at a mean z — 0.97. Using wider projected radius bins we find that 
for the galaxies with redshifts 1.2—1.6, £(r) ~ (r /1.0e ±0 2 h~ l Mpc) -1 ' 8 at a mean z ~ 1.39. This 
high redshift correlation is formally an upper limit on the basis of the bootstrap error estimates, 
although the Poisson errors do indicate a significant measurement. We will accept the result as an 
indicative measurement, because it remains unclear how best to estimate the errors. The galaxies 
in this highest redshift bin are dominated by the faint B subsample and are not complete in the 
K subsample This selection bias likely leads to an underestimate of the correlation at this redshift 
if the enhanced correlation of red-selected galaxies over blue-selected ones seen at lower redshift is 
present at this redshift. 

The decrease in the correlation length with increasing redshift is significant within these data 
although the quoted errors are purely the internal errors of the fit and do not account for field-to- 
field differences beyond these two patches. The fitted e = 0.2 ± 0.5, based on our K-hand sample 
alone. 



5. Correlation Dependence on Color and Luminosity 



The color dependence of faint galaxy clustering is an important clue to the formation mech- 
anisms of galaxies. Higher mass galaxies are expected to be more correlated at formation and 
subsequently environmental effects can modify the correlations {e.g. Loveday et al. 1995| ). In 
Figure || the w p {r p ) are shown for the 0.3 < z < 0.9 galaxies with colors redder or bluer than 
(U — K)q = 2 mag (k-corrected, rest-frame ratios between the flux f v at 3500A and at 21000A on 
the AB magnitude system, see Cowie et al. 1996). The red galaxies, with ro = 2.4e ±0 ' 14 /i~ 1 Mpc, 
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are about 5 times (with about a 50% error) more strongly clustered than the blue objects, which 
have ro = 0.9e ±0 ' 22 fo _1 Mpc, both at fixed 7 = 1.8. The red galaxies appear to have a substantially 
steeper correlation slope, 7, than our adopted value of 1.8, although a survey with larger sky area 
is needed to assess the slope. The average Mk = —22.5 mag for the blue sample, whereas it is 
Mk = —23.9 mag for the red sample, which is quite a small luminosity difference for the large 
correlation difference ( Loveday et al. 1995 , Tucker et al. 1996| ) . The blue population, at a mean 



B of 23.7 mag, has strong [0 11] lines, with a mean W eq of 41A. Therefore, these are very likely to 
be substantially brightened relative to their intrinsic luminosities at a more normal star formation 
rate. The blue galaxy population must largely disappear from the normal galaxy population at 
low redshift, since its correlation length would grow to only about 2 h -1 Mpc given our estimated 
e ~ 0.2. Within the luminosity range explored in the low redshift APM survey, no population is 



this weakly correlated (Loveday et al. 1995) although the low luminosity, M& = —18.6, late type 
galaxies have tq = 2.9±0.4/i _1 Mpc. We conclude that these high redshift faint blue galaxies cannot 
become a significant component of the low redshift galaxy population above luminosities of 0.4L*. 
It is possible that these objects merge with higher luminosity galaxies to drive their evolution. 

The clustering of lower luminosity galaxies is weaker than that of high luminosity galaxies. 
For 0.2 < z < 0.4 the galaxies with M K > -21.5 have r = 1.8e ±a2 /i _1 Mpc and for 0.4 < z < 0.8 
those with Mk > —23.5 have ro = l.le ±0 ' 2 /i _1 Mpc. In both these redshift ranges, the higher- 
luminosity galaxies are more correlated that the lower luminosity ones. The luminosity difference is 
about 2.5 magnitudes. It is quite unlikely that the high-redshift, high-luminosity, weakly-correlated 
population could have evolved into a weakly-clustered, low-luminosity, low-redshift population. 
First, it requires that the clustering is being reduced in physical co-ordinates with time. Second, 
about 4 magnitudes of fading per unit redshift is required, whereas a much slower rate of luminosity 
density evolution of the population average is observed ( Lilly et al. 1995| ). Third, it would require 



that the high-luminosity galaxies at lower redshift originate from some unobserved high redshift 
population — an effect which ET-band observations minimize through their relatively small redshift 
corrections. On the other hand, this discussion does emphasize that a more precise definition of 
comparable galaxy samples at low and high redshift would be very desirable, ideally done based on 
the intrinsic properties of the galaxies themselves. 

The cross-correlation of low-luminosity galaxies, Mk > —23.5 mag, with high-luminosity galax- 
ies, Mk < —23.5 mag, is shown in Figure [?] for the 0.3 < z < 0.9 range (expanded slightly to boost 
the sample size). Although the sample is smaller than is really desirable it shows quite intrigu- 
ing how much more strongly the low luminosity galaxies cluster to high luminosity "hosts" within 
lOO/i -1 kpc, beyond which the cross-correlation drops to a value similar to the field cross-correlation 
of low- luminosity galaxies. The effect is directly visible in the redshift diagrams of Figure ||| The 
enhanced cross-correlation of close pairs is not primarily the result of an enhanced star forma- 
tion which would raise the optical band luminosity between 1 and 2 magnitudes, and hence would 
increase the numbers above some flux limit. The lower luminosity galaxies have substantially 
stronger [On]; however the expected accompanying increase in the K luminosity is small. Conse- 
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quently these data favour a genuine increase in close pairs at small separations over the power-law 
correlation. 



6. Evolution of Galaxy Correlations 

To follow the clustering amplitude as a function of redshift we need to compare the results of 
a variety of redshift surveys, having varying r$ and 7. The quantity Tq , is a useful measure of the 
clustering amplitude, which can be interpreted either as the amplitude at 1 /i -1 Mpc, or, given that 
these quantities are normally the results of a fit to data over a range of scales, is a measure of 
the average correlation within a fixed proper volume. 

At low redshift there are several recent measurements of clustering with a range of sample 
definitions. The 6j-selected APM survey finds tq = 5.1 ± 0.2/i _1 Mpc and 7 = 1.71 ± 0.05 giving 



16.2±2.5 at z ~ 0.06 QLoveday et al. 1995|) . The ^-selected LCRS finds r = S.OiO.U/i" 1 Mpc 



and 7 = 1.79 ± 0.04 giving rj = 17.8 ± 2.2 at a mean z ~ 0.1 ([Lin 1996) ) . The IRAS-selected 



correlation function is tq = 3.76 ± 0.2/i _1 Mpc and 7 = 1.66 ± 0.11 ( [Fisher et al. 1994j ) for 
Tq = 9.0 ± 2.6. There is no large -KT-band selected redshift survey at low redshift; however the 
angular correlation of bright K-selected galaxies is measured, from which we estimate Tq = 27.5 
and 13.8 at z = 0.13 and 0.23, respectively, with ~20% systematic errors ([Baugh et al. 1996| 



referred to as BFFS). A preliminary measurement in the LCRS survey finds that the red galaxies 
have a correlation comparable to these /f-band results ( Tucker 1994, 1996| ) . 



At higher redshift we have the results of this paper, the r-selected CNOC field survey at a mean 
redshift of 0.36 ( [Shepherd et al. 1997 and work in progress), which gives Tq ~ 9.6 for the red half of 



the sample and 7.5 for the blue half. The CFRS results flLeFevre et al. 1996| ) provide correlation 
estimates over the redshift range 0.2-0.9. They find £(r) = (r/2.4±0.17/i _1 Mpc) -1 ' 64 for 0.2 < z < 
0.5, £(r) = (r/1.4 ± O.Wh' 1 Mpc) -1 ' 64 for 0.5 < z < 0.75, and f (r) = (r/1.4 ± 0.20/1" 1 Mpc)" 1 - 64 
for 0.75 < z < 1, where we have adjusted their correlation lengths for qo = 0.5 to our qo = 0.1 using 
their formula. The CFRS data indicate 7 = 1.64, which would be a relatively poor description 
of the -RT-selected data here. As a direct comparison with the CFRS measurements, we define a 
complete subsample limited at I = 22.5 mag from our data. We find that r = 1.9e ±019 /i _1 Mpc 
for 0.3 < z < 0.9, which is statistically identical to the CFRS measurement of 1.8 ± 0.18/i _1 Mpc 
(adjusted to qo = 0.1) over the < z < 1 range. 

The correlation amplitudes are plotted against redshift in Figure || For comparison the simi- 
larly measured correlations from n-body simulations ( |Colin, Carlbcrg fc Couchman 1996j ) are also 
shown, all scaled with a linear multiplicative factor roughly to fit the LCRS correlation measure- 
ment. The fitted 7 values from these simulations are in the range 1.8-2.0, which is compatible with 
the .fT-band sample but not the shallower slopes usually seen in optically selected samples. The 
scaling factors are in the range of 0.55 to 0.70, which is a little larger than desirable (ideally one 
would do specially matched simulations), however, these factors are small compared to the factor 
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of 20 or so in the evolution of the correlation functions. 

There are two conclusions to be drawn from Figure First, optically selected galaxies appear 
to always be significantly less correlated than iT-selected galaxies, the difference being typically 
about a factor of two in the amplitude. Second, for the iT-band and red-selected samples, which 
should be least sensitive to galaxy population evolution, we see that the evolution of galaxy clus- 
tering is reasonably well described by an Q = 0.2 model. However, the amplitude measured in the 
n-body simulation is multiplied by a factor of 0.55, which is approximately the square of the bias 
factor, b = 0.75. Hence, either the galaxies are anti-biased with respect to the matter clustering, or, 
the normalization of used for the n-body simulations should have been approximately <7g — 0.75. 
The observed i^-band selected evolution appears to be too slow with redshift to readily agree with 
the = 1 predictions, although the formal level of exclusion is strongly dependent on the low red- 
shift normalization. The second conclusion to be drawn is that the amplitude of the correlations 
increases as the color used to select the objects becomes redder and there is weak evidence that 7 
is also steeper in the K sample. 



7. Conclusions 

The fitted correlation length of luminous X-selected galaxies over the redshift range 0.2 to 
1.2 is substantially stronger than that found for optically selected samples, about a factor of two 
in the amplitude, r$. The galaxy correlation amplitude is measured at a mean z ~ 1.39 as tq = 
1.0e ±0 - 2 /i _1 Mpc (formally an upper limit, but deficient in the more strongly clustered faint red 
galaxies), r = 1.4e ±0 - 15 /i _1 Mpc at z ~ 0.97, 2.0e ±0S5 h~ 1 Mpc at z ~ 0.62, and 2.9e ±0S2 h~ 1 Mpc 
at z ~ 0.34. Together these give a clustering e ~ 0.2 ± 0.5. 

The red galaxies are about a factor of 5 more correlated than the blue galaxies, which have 
ro — 0.9e ±0 ' 22 /i -1 Mpc. These blue galaxies have a mean equivalent width in the [On] line of 4lA. 
Together this can be taken as strong evidence that the faint blue galaxies are an intrinsically weakly 
correlated population (therefore likely low mass) with a high star formation rate that brightens 
them into the range of much more strongly correlated red galaxies. These galaxies are so weakly 
correlated that for our measured growth of correlations, e ~ 0.2 ±0.5, they would grow to a current 
epoch correlation length of ro — 2h~ l Mpc. This correlation length is less than that measured in 
any galaxy population at low redshift ( Loveday et al. 1995| ), so these faint blue galaxies cannot by 



themselves make a significant contribution to the current epoch galaxy population. 

Overall the K selected galaxy correlation evolution is somewhat too slow with redshift to be 
easily consistent with the evolution of the matter correlation function for f2o — 1- Affs- 0.8 and 
Q ~ 0.2 — 0.3 would describe both the amplitude and its evolution, if these galaxies are tracing the 
matter clustering. To further test models of correlation evolution requires large datasets of precision 
comparable to that available in current low redshift surveys with good control over population 
changes with redshift. 
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Fig. 1. — The combined redshift distribution of the 248 galaxies with measured redshifts in the 
B,I,K magnitude-selected sample of the SSA13 and SSA22 fields. The bins have Az = 0.01. 

Fig. 2. — The positions of if -detected objects on the sky in SSA13 and SSA22. The small ticks 
are at intervals of 20" in RA and 10" in Dec. Symbol area is proportional to txik- Galaxies are 
shown as squares, stars as crosses, objects without confident redshift identifications as triangles, 
and unobserved objects as plus signs. A larger fraction of objects in the SSA22 field are stars due 
to its lower galactic latitude. 

Fig. 3. — The redshift vs projected proper distance from field center in the RA direction, calculated 
for qo = 0.1 for the entire sample. The area of each circle is proportional to the object's if-band 
luminosity. The size for a galaxy with Mk = —23.5 is shown near the bottom of the plots at 
z = 1.5. 

Fig. 4. — The angular correlation of the photometric sample (including stars). The two points at 
largest separation have negative values, which are not statistically significant. The indicated la 
errors are from a bootstrap analysis with 100 resamplings. 

Fig. 5. — The projected real space correlation function for redshift subsamples at 0.2 < z < 0.4 
(top-left panel), 0.4 < z < 0.8 (top-right panel), 0.8 < z < 1.2 (right panel), and 1.2 < z < 1.6. 
The la bootstap errors (narrow error flags) indicate that the high redshift correlation is an upper 
limit, although the Poisson error bars (wide error flags) indicate a significant result. The high 
redshift subsample is deficient in faint red galaxies, which are expected to be strongly correlated. 
The point at 0.01/i _1 Mpc corresponds to an angle of about 2" where galaxy images overlap and 
is not used to calculate the correlation length. All errors are from a bootstrap analysis and are 
considerably larger than Poisson estimates. 

Fig. 6. — The projected real space correlation function for the blue and red subsamples over 
0.3 < z < 0.9. Bootstrap errors are shown. 

Fig. 7. — The projected real space cross-correlation function of the Mk > —23.5 subsample with 
the Mk < —23.5 subsample. The line shows the auto-correlation of high luminosity galaxies over 
the same redshift range. Bootstrap errors are shown. 

Fig. 8. — The evolution of the nonlinear amplitude of the correlation function, Tq. The data are 
compared with results from other surveys in the literature (see Section 6 for details), with all 
measurements adjusted to go = 0.1. The results from this paper are shown as filled circles, with 
superposed crosses for the blue and red subsamples. The solid, dashed, and dotted lines are based 
on power law fits to the £(r\z) measured in n-body simulations of CDM universes for different 
assumed values of and A, renormalized to pass through the LCRS measurement at low redshift. 
Note that correlations of relatively red selected galaxies are always larger than those of blue selected 
galaxies. Bootstrap errors are shown. The Keck and BGFS samples are K selected. CNOC and 
LCRS are r selected, and CFRS is / selected. The APM is blue selected and the IRAS sample is 
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dominated by relatively blue galaxies which contain warm dust. 
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